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ABSTRACT 

There are many fluid flows where the onset of transition can be caused by different 
instability mechanisms which compete in the nonlinear regime. Here the interaction of 
a centrifugal instability mechanism with the viscous mechanism which causes Tollmien- 
Schlichting waves is discussed. The interaction between these modes can be strong enough 
to drive the mean state; here the interaction is investigated in the context of curved 
channel flows so as to avoid difficulties associated with boundary layer growth. Essentially 
it is found that the mean state adjusts itself so that any modes present are neutrally 
stable even at finite amplitude. In the first instance, the mean state driven by a vortex 
of short wavelength in the absence of a Tollmien-Schlichting wave is considered. It is 
shown that for a given channel curvature and vortex wavelength there is an upper limit 
to the mass flow rate which the channel can support as the pressure gradient is increased. 
When Tollmien-Schlichting waves are present then the nonlinear differential equation to 
determine the mean state is modified. At sufficiently high Tollmien-Schlichting amplitudes 
it is found that the vortex flows are destroyed, but there is a range of amplitudes where 
a fully nonlinear mixed vortex-wave state exists and indeed drives a mean state having 
little similarity with the flow which occurs without the instability modes. The vortex and 
Tollmien-Schlichting wave structure in the nonlinear regime has viscous wall layers and 
internal shear layers; the thickness of the internal layers is found to be a function of the 
Tollmien-Schlichting wave amplitude. 
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§1. Introduction. 

Our concern is with the strongly nonlinear interaction of large amplitude Tollmien- 
Schlichting waves and Gortler vortices in curved channel flows. Interaction problems of 
this type are relevant to many flows of practical importance where viscosity and centrifugal 
effects combine to stimulate travelling wave disturbances and stationary vortices. The most 
obvious examples are the flow driven by a steady pressure gradient in a curved duct and the 
flow over a laminar flow aerofoil, see Harvey &: Pride (l). In the latter problem, a question 
of fundamental importance is whether the presence of finite amplitude Gortler vortices and 
Tollmien-Schlichting waves and their consequent interaction cause the premature onset of 
transition. 

Here, because of the difficulties associated with boundary layer growth, we concentrate 
on fully developed curved channel flows driven by azimuthal pressure gradients. We shall, 
however, indicate the probable relevance of our calculations to the external flow situation. 
The interaction which we consider is strong enough for the Tollmien-Schlichting waves 
(hereafter referred to as TS waves) to have an 0(1) effect on the steady vortex flow which 
exists in the absence of the waves. Likewise this 0(1) correction to the vortex state has 
an 0(1) effect back on the wave so that the vortex and the TS states are strongly coupled. 
Remarkably this strong interaction occurs at extremely small TS amplitudes; it is the large 
inertial effects associated with both the TS wave and Gortler vortex which generate this 
large response. 

The basis for the expansion procedure to describe the interaction is that given by 
Hall & Lakin (2) for small wavelength Gortler vortices in growing boundary layers. The 
nonlinear states considered in that paper have the ‘mean’ state driven by the vortices 
which are themselves trapped between viscous shear layers in which the vortex amplitude 
decays to zero as a solution of a nonlinear Airy equation. We shall apply this theory to 
channel flows but with the complication that the mean state is simultaneously driven by 
the vortices and the TS waves. The ‘core equation’ to determine the mean state has a 
significantly different structure in the presence of a TS wave. An important consequence 
of this change in structure is that vortex states cannot exist beyond a critical TS wave 
amplitude. Before discussing our problem in more detail, we shall give a brief review of 
related work on TS-Gortler vortex interactions. 

Most of the previous work in this area has in fact been in the context of fully-developed 
flow. However, one of the earliest papers, Nayfeh (3), did consider the external flow case 
but there the amplitude of a linear Gortler vortex was simply imposed and the linear insta- 
bility of the resulting flow to an infinitesimal TS wave was considered. The linear Gortler 
vortex was obtained by ignoring boundary layer growth; since the vortex wavelength cho- 
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sen was comparable with the boundary layer growth rate we know from the work of Hall 
(4) that the eigenfunction used by Nayfeh does not satisfy the linear stability equations. 
Furthermore, the amplitude of a Gortler vortex must be calculated from the Navier-Stokes 
equations using the approach of Hall (5), Hall (6) or Hall & Lakin (2). Having made the 
simplification discussed above, Nayfeh found that small vortices have a large effect on TS 
growth rates. Later Malik (7) repeated Nayfeh’s calculations and showed that Nayfeh’s 
conclusions were incorrect and indicated a numerical error in Nayfeh’s work. 

Subsequently, Bennett & Hall (8) investigated the linear instability of nonlinear vortex 
flows in curved channels to both 2D and 3D TS waves. Here the nonlinear states were 
found by solving the full Navier-Stokes equations numerically. The particular TS waves 
imposed on this flow correspond to the lower branch of the neutral curve. This type of TS 
wave is particularly relevant in external flows, Smith (9,10), so it is of course important 
to understand how vortex flows influence this mode even though more unstable modes 
at finite Reynolds numbers can occur in channel flows. Bennett & Hall (8) showed that 
vortex flows with wavelengths comparable with the channel width can have a significant 
destabilizing effect on lower branch TS waves. In a subsequent paper Bennett, Hall & 
Smith (11) extended the work of (8) to the nonlinear stage where the TS wave has an 
O(l) effect back on the vortex flow. The vortex wavelengths are then comparable with 
the channel width and the flow can only be determined by numerical calculations. The 
somewhat limited number of calculations performed in (11) showed that nonlinear effects 
prevent the exponential growth of the TS amplitude and that a supercritical bifurcation 
to a mixed TS-vortex state occurs. 

If the vortex wavelength is taken to be long (in fact comparable with the wavelength 
of the TS wave), a different type of interaction takes place. In this case, Hall &: Smith 
(12) show that initially the interaction is governed by amplitude equations appropriate 
to a resonant triad. The equations were found to possess a finite time singularity whose 
structure indicated a more nonlinear state. The latter state is governed by a linear partial 
differential equation for the vortex coupled to an integro-differential amplitude equation 
for the TS wave. It was again found, even in the more nonlinear state, that a finite 
time singularity exists. In a related paper for channel flows, Hall & Smith (13) showed 
that the vortex velocity field develops a shorter spanwise scale where certain types of 
singularity occur. Thus, it is important to understand the vortex-TS interaction at smaller 
wavelengths comparable to or smaller than the channel width since there is some evidence 
from the work of Hall & Smith (13) and some related boundary layer calculations of Hall 
& Smith (14) that these shorter scales are induced by an interaction involving longer 
wavelength vortices. 
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Finally, in our discussion of previous work we mention the calculation of Daudpota, 
Hall & Zang (15). This interaction problem concerns small amplitude TS waves and 
vortex states at finite Reynolds numbers at a value of the curvature at which both modes 
are equally likely. The results of that paper suggest that Gortler vortices can have a 
significant effect on the equilibration of TS waves. 

In the present paper, we consider the interaction of small wavelength Gortler vortices 
and lower branch TS waves. Our reason for concentrating on these limits is easily explained 
by reference to the corresponding external flow case. There it is well known that Gortler 
vortices set up in an experiment conserve their wavelength as they move downstream. 
Since the boundary layer itself thickens it follows that the local nondimensional vortex 
wavenumber becomes large as the vortex develops downstream. Thus, the small wavelength 
limit in the external Gortler problem is appropriate to the ultimate development of any 
initial fixed wavelength vortex. Similarly, we concentrate on lower branch TS waves since 
they are the most unstable modes which can be described by formal asymptotic methods. 

Since there has been no previous work on small wavelength large amplitude vortices in 
channel flows, we extend the work of Hall & Lakin (2) to fully developed flows. We use the 
condition of constant mass flow rather than constant pressure gradient and obtain curves 
of mass flow against pressure gradient in the presence of vortices for different values of the 
curvature. We show that when vortices develop they cause unusually large increases in the 
pressure gradient to be induced in order to produce a small increase in the flow rate. 

We then consider the effect of a small amplitude TS wave on the above flow structure. 
Initially the vortex wavelength is taken to be comparable with the channel width so that 
the interaction equations of Bennett et al. (11) apply. By taking the further limit of 
small vortex wavelength, we derive an asymptotic description of a fully nonlinear flow 
with ‘mean’ part simultaneously driven by TS waves and Gortler vortices. 

We show that vortices can only exist in a finite region bounded by thin shear layers 
whose structure depends on the size of the TS wave amplitude. We will see that supercrit- 
ical states exist in the presence of TS waves and that, at a given Taylor number, several 
equilibrium states with different TS amplitudes and frequencies are possible. We shall 
show that there is an upper bound for the TS wave size at a given Taylor number. The 
procedure adopted in the rest of this paper is as follows: in §2 we formulate the nonlinear 
equations governing 3D flows in curved channels. In §3 we solve these equations for strongly 
nonlinear vortex flows and in §4 we extend this to allow for the simultaneous presence of 
TS waves. Finally, in §5 we discuss the results of §3,4 and draw some conclusions. 
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§2. Formulation of the problem. 


Consider the flow of a fluid of density p and kinematic viscosity 1 / in a curved channel 
with walls defined by r = Ri and r — R 2 (> Ri) with respect to the usual cylindrical 
polar co-ordinates (/, 0' , z). We assume that the curvature of the channel 6 = (R 2 — Ri) / R\ 
is small. The basic steady flow considered is that of an azimuthal velocity field driven by 
an appropriate pressure gradient. It is convenient to define variables x, y, z, t by 

RiO' _r-R 1 _ z , _ U 0 t' 

X R 2 — R\ V R 2 -R^ Z R 2 — R\' Re(R 2 - Rj ’ [ ’ 

where Uq is a typical azimuthal velocity, t is a dimensionless time and Re = Uod/v is 
the Reynolds number. Writing the velocity field scaled on Uq with respect to (x, y, z) as 
(u,v,w), scaling the pressure on pU^ , and substituting into the continuity and Navier- 
Stokes equations yields the system 


1 du Su dv dw 

H "TT + ^ 1" "5 — = 

F dx F dy dz 


1 

Re 

1 

Re 


dt F 1 


26 dv 1 dp 6uv 

u + — — — — — — = Nu H , 

ReF 2 dx F dx F 


V 2 - - 

dt 

1_ 

Re 


6 2 _ 
V 2 


v — 


26 du 

rJf 2 ^ 


dp_ M _Su 2 
dy F 


d_ 

dt 


dp 

w — — = Nw, 

dz 


where 


F = 1 + 6y, 

and the operators V 2 and N are given by 

1 


V 2 = 


6 d 


d 2 d 2 d 2 
F 2 d^ + dy 2 + dz 2 + Fdy ’ 


.. u d d d 

N = 1- v b w — . 

Fdx dy dz 


(2.2a) 

(2.26) 
(2.2c) 
(2. 2d) 

(2.3a) 

(2.36) 

(2.3c) 


The Taylor number T is given by T = 2 Re 2 6, and here we take this to be an 0(1) 
parameter. Dean (16) demonstrated that instability of the basic flow given by 


p=-^ + 0(6), u= y(y-y 2 ) + 0(6), (2.4) 

(where D is a constant) occurs first for 0(1) axial wavelengths and at an 0(1) Taylor 
number so that here the Reynolds number is large. 
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We now consider the case of fully nonlinear small wavelength Gortler vortices in the 
channel and suppose that the mean flow correction to the basic flow in the presence of the 
vortices is as large as the basic flow itself. 

§3. The case of a fully nonlinear Gortler vortex in the channel. 

Here we follow the ideas of Hall & Lakin (2) who studied the problem of fully nonlinear 
vortices in an external growing boundary layer. The necessary scalings for the present 
problem follow those given in that paper with a related flow structure. This structure 
consists of a vortex activity region (7) defined by (0 <)yi < y < yi{< 1) say, with thin 
layers ( 77a, 6 ) surrounding y = yi and y = y 2 which act to smooth out the algebraically 
decaying vortices in the central region. The structure is completed by regions ( IIIa,b ) 
defined by 0 < y < yi and < y < 1 in which, to leading order, the flow is merely driven 
by a constant pressure gradient with no vortex present. 

Unlike Hall & Lakin (2), here our flow velocities (u,v,iu) are independent of x and so 
we seek steady solutions which assume the forms 

[«,v,w) = (u(y,z),± ; V(y,z),± ; Wh,z)), P=-ff+^^> (3-1) 

in the notation of (2.2). Substituting (3.1) into the system (2.2) and ignoring small terms 
yields, at leading orders, the governing equations 


dV 

aw 


(3.2a) 

dy + 

dz 

= 0, 

b 

i 

i 


a 2 u a 2 u 

(3.26) 

V dy +W dz 

■ D+ ay’ + *.■ - 

dV dV TU 2 


ap ( a 2 v t a 2 v 

(3.2c) 

v a y +w a, + 2 

' — 

dy dy 2 * dz 2 ’ 

T aw aw 

dP 

d 2 W d 2 W 

(3. 2d) 

V 6y +W e,=- 


dy 2 ' dz 2 


We look for a vortex with small axial wavelength /? C 1, and use the results of (2) that 
the Taylor number T expands in the manner 

T = To/? -4 + Tx/T 3 + • • • , (3.3) 

and that in the main vortex activity region (7) the appropriate velocity and pressure 
expansions take the forms 


U = u 0 + pu l + ... + {/3 E(Uq + pUl + ...)+ p 2 E 2 {U 2 + ...) + ... + c.c}, (3.4a) 


5 


V = {p- l E{V£ + pvt + ...)+ E 2 (Vo 2 + + C.c}, (3.4 6) 

W = {E{w£ + (3Wl + ...) + 0E 2 (W* + ...) + ... + c.c}, (3.4c) 

P = P~*Po + P'*Pi + • • • + {/T^Po 1 + ppl + ...)+ P 2 (P 2 + ...) + ...+ c.c}, (3.4d) 

where E = exp(iz/ (3) and c.c denotes the complex conjugate. Notice that in the absence 
of the vortex the flow reduces to the steady state with V = W = P = 0 and U = y(y — y 2 ). 

The expansions (3.4) are substituted into (3.2) and powers of /? are equated for each 
Fourier coefficient. The zeroth order equations for the vortex in the core are 


HV 1 

W +iWi = 0 ’ 

TJ 1 + v -1 — - 0 
u ° + 0 dy ~ ’ 

(3.5a, 6) 

FJ) 1 + T 0 UqU 0 — 0, 

t'Po 1 + Wl = 0, 

{3.5c, d) 

and consistency of (3.56, c) flow in the core 



dUn 


(3.6a) 

ToUo-r?- 

dy 

= 1. 


Further, the vortex amplitude | V 0 X | follows from the x- momentum equation (3.26) and 
satisfies 


d?Un d 

0 = -D- 2 — 




(3.66) 


dy 2 ~ ~dy [dy 

Consequently, the equation (3.6a) for the mean flow fixes the vortex. Integrating (3.6a) 
gives 

Bo = ,/fe+5l. (3.6c) 


To 

where a is a constant and then (3.66) yields the condition 


du n 

-y— = 6 — Dy 
dy 


1 2 

2 dy |K ° 1 ’ 


(3.6d) 


where 6 is another constant. As IVqI 2 is non-negative, the extent of the vortex activity 
region is defined by | VJj 1 1 = 0. If y = yj and y = yi(> yi) are these edges of activity then 

1 


6 = Dyj + 


\/T 0 (a + 2 y,-) 


(3.7) 


for j — 1,2. The thicknesses of the thin shear layers centered on yi and y 2 are found by 
balancing diffusion across the layers with convection in the streamwise direction and are 
found to be 0(/?t). Hence, in the upper layer, we define the 0(1) co-ordinate f by 


e = 


y - y 2 

/?! 


(3.8) 
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The expansions in this region ( lib ) are then 

U = Uq + + /?U2 + P*u$ + . . . + { 0 *E(Uoi + f}*Uu + . . .) P*E 2 (Uq2 + ...)+...+ c.c}, 


( 3 . 9 a) 

V = {p~lE{V 01 + piv n + ...)+ P*E 2 {Voz + ...) + ••• + ™}, ( 3 . 96 ) 

W = {P~sE{W 01 + 0 *Wu + ...) + PE 2 (W 02 + ...) + ... + c.c}, ( 3.9 c) 

P = P P 8 Pj + . .. + {/? *^(Pqi + /? 8 -Pn + • • •) +P^E 2 {Pq2 + • • •) + • • • + c.c}. ( 3 . 9 d) 


The governing equations are very similar to those presented in ( 2 ) and so we will merely 
summarize the important features of their solutions. The leading order terms in the mean 
flow expansion satisfy 

d 2 U o d 2 U i 

He ~ 1? ~ ’ 

and matching with the vortex activity region yields 


u 0 


-ft- 


+ 2y 2 ) 


Ui = 


^/T 0 (a + 2y 2 ) 


The governing equations for (£/oi, V01) are found to be always consistent, but those for 
Un and Vn admit a solution only if a certain solvability requirement is met. This condition 
is obtained on consideration of higher order mean flow terms and leads to a problem (which 
is a particular form of the second Painleve transcendent) given by 


d 2 V 01 

d? 


+ £ 


[ 3 (a + 23/2) 3 


— — \/To(a + 2y 2 ) 


2 „ \/ro(a + 2y 2 ) 

Vbi = \V 01 \V Ql \ 2 - V ° V 3 — /V 01 , 


where / is some constant whose actual value is not required here. This equation has been 
shown ( 17 ) to have a solution for which 
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2 


D\Jr^a + 23 / 2 ) — 


as £ — ► —00, 


|F 0 i| 2 — ► 0 as £ — ► 00. 


Thus the fundamental terms Uq\ and V01 decay exponentially to zero as £ — ► 00 and 
so the vortex is constrained to lie below 3/2. A completely similar analysis holds for the 
lower shear layer (I I a) at y = 3/1. We see that the mean flow terms play an essentially 
passive role in these shear regions, in as much as that the leading order mean flow term 
and its derivative must be continuous across the shear layers. Outside the vortex activity 
zone in regions (Ilia, b) we have 


t/ = § + 0(/^), 
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and (3.66) simplifies to the governing equation = — D. 

Summarizing at this point then, we now have (given To) sufficient information so as 
to enable us to determine the shear layer positions j/ x and y 2 and the constants a, 6 in 
(3.6c, d). The leading order mean flow quantity u is governed by the following constraints 
and conditions: 

Zone (I). 

' ( “ + 2!,) (3.10a) 


u = 


To 


Zones (Ilia, 6). 
Here 


d 2 u 


(3.106) 


and u and ^ are continuous across the shear layers (Ila, 6) at y\ and y 2 . From (3.7) we 
also have 

Dyi H — , = Dy 2 H — , , (3.11) 

y To(a + 2 y x ) y To(a + 2y 2 ) 

and no-slip conditions u = 0 on y = 0, 1 are imposed. 

The constant D is not determined as yet and we choose it so that the non-dimensional 
mass flux across the channel is kept constant for varying values of To. We take the nor- 
malization such that u = y — y 2 in the absence of any vortices so that D is fixed so 

1 


/ 


udy = 


6 


(3.12) 


We used a simple iterative technique to solve the problem (3.10), (3.11) and it was 
found that the method was most easily formulated in terms of the quantity TqD 2 . Given 
this value, a guess was made for the shear of the mean flow at the lower wall y = 0 and 
(3.106) was integrated to the position y\ where To«^ = 1. Continuity of the mean flow 
across the shear region (Ila) yielded a value of a on using (3.10a), y 2 was then determined 
from (3.11) and finally (3.106) integrated to y = 1. If u(l) did not vanish, the assumed 
shear at the lower wall was adjusted until this condition was met. Finally, D was fixed so 
as to satisfy (3.12). 

The results of these calculations are shown in figures (l) & (2). The work of Hall (5) 
demonstrated that on the basis of linear theory a vortex can exist for T 0 = T c = 6\/3. 
Not surprisingly, our computations revealed that a solution of the problem (3.10)-(3.12) is 
only possible for T 0 > T e . Also, no solutions occur for T 0 > 32 and, indeed, as T 0 — * 32, 
the mass constraint (3.12) implies that the pressure gradient D — ► oo. Figure (l) shows 
the positions of the shear layers yi, y 2 as a function of T 0 and as T 0 — * 32 these positions 
migrate to the edges of the channel, whereas as T 0 tends to the linear value T c they merge 


8 



at y = |(1 — ^-). Figure (2) illustrates the dependence of the Taylor number T 0 on the 
pressure gradient D. 

The normalization we have used to determine D in our calculations thus far is that to 
ensure a constant non-dimensional mass flow rate. In a realistic experiment, it would be of 
more interest to determine the dimensional mass flow rate as a function of the total pressure 
drop along the channel. On recalling the non-dimensionalizations introduced previously, 
we may easily translate our results so that they may be interpreted in dimensional terms. 
Figure (3) illustrates the dependence of the dimensional mass flow rate on the pressure 
drop along the channel. In the absence of any vortices, we know that the mass flux is 
proportional to the pressure drop, and it is seen that as vortices develop, a large increase 
in pressure drop is required to produce a small increase in the flow rate. The precise broken 
curve on figure (3) which is actually followed in any particular experiment is determined by 
the physical parameters of that experiment, and we see that in each particular case there 
is a maximum flow rate that can be achieved in the presence of the vortices, irrespective 
of the size of the pressure gradient. 

Lastly in this section, we see that two asymptotic regimes naturally arise here. Firstly, 
forT 0 D 2 = 24v/3 + £,0 < 5* -C 1, i.e. returning to the case of very small vortices, we find 
that 



and To = 6\/3 + (S*/2) + — These results imply the scalings required for a weakly 
nonlinear analysis of the vortices along the lines of Hall (18). 

Secondly, for 0 = TqD 2 1, we obtain 


1 24 

Vl — T + TT + • • • > 

0 04 




21 


a a 

0 04 

To = 32 - 2^0-7 + 2*0"* + . . , 


This demonstrates that To < 32 in all cases, a restriction which is a direct consequence 
of the constant non-dimensional mass flow constraint. As T 0 approaches this upper bound 
the vortices occupy nearly all of the channel and the mean flow profile takes a square root 
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form everywhere except for thin regions at the sides of the channel. If the non-dimensional 
flow constraint is relaxed and replaced by a constant pressure gradient requirement, the 
vortex may then persist for arbitrarily large values of To. 

The effect on the vortices of introducing TS waves into the channel is now investigated 
and the scalings are chosen so as to ensure a strong interaction between the vortices and 
the waves. 


§4. Strong interactions between small wavelength Gortler vortices and lower 
branch TS waves. 


We now extend the analysis presented in the last section to the case when TS waves of 
azimuthal wavenumber O (Re~$) are also present in the flow. The size of the TS wave is 
chosen to be such that the nonlinear terms in (2.2) arising from interactions of the TS wave 
with itself, are of comparable size with the leading order vortex terms. These interactions 
are most important in the core of the flow away from the viscous wall layers in which the 
TS wave adjusts to the no slip conditions required at the walls. 

Bennett et al. (11) examined the problem under consideration here but took the vortex 
and TS wave to vary on an 0(1) scale in the 2 -direction, whereas we choose to study the 
case of a small axial lengthscale. The TS wave is assumed to be proportional to 

E = exp |t — J dt^j j , (4.1) 

where the wavenumber a and the slowly varying frequency fi are real, and where, antici- 
pating the ensuing analysis, it is found convenient to define the small parameter e = Re~$. 
Away from the viscous wall layers at y = 0 and y = 1 we follow (11) and consider velocity 
and pressure fields which take the forms 


u = u 0 + (e e uiE + c.cj (l + 0(e 7 )j , 

(4.2a) 

v = 2 + (e 7 t/i E + c.c) (l + 0(e 7 )) , 

(4-26) 

w = + (e 7 w iE + c.cj (l + 0(e 7 )) , 

(4.2c) 

—Dxe 7 + + (e 8 PiE + c.cj (l + 0(e 7 )j , 

(4.24) 


where the unknown quantities are independent of x. 

Substitution of these forms into the governing equations (2.2) yield that at leading 
orders the vortex and TS parts of (4.2) satisfy 



dw o 

~d7 


= 0 , 


(4.3a) 
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(4.36) 


and 


( 


(— 

\dy* 


d 2 

d 2 

d\ 

| uo + D - 

1 

( du o 


du 0 \ 

dy 2 + 

dz 2 

dt) 

” 2 

V°dy 

+ W 0 

~d7j 

d 2 

d' 

l 

T 2 

dpo 

if 

dv a 

+ U>o 

+ d^ 

dt j 

J v 0 

2 U ° 

dy 

= 2 r * 

+2 


► duj 

dv:\ 

+ u, ‘ar) 

+ c.c 

9 


dv { 

~dz 


r) 


(&_ d*_ 

\dy 2 dz 2 


d \ dp 0 1 ( dw 0 dw 0 \ 

m) w °--d 7 = 2r a7 + u;o 17j 


+2 


K dw\ dw 

- tau 1 w 1 + Vx~^ + Wi 


dy 


hv\\ 
dz ) 


+ C.c 


dvi dwi 

tau l + 1 - “5 — — 0, 

ay dz 

du 0 du 0 

iau 0 Ui + — Vi + — u>i = 0, 


tau 0 Vi = 


iau 0 Wi = 


dp i 
dy’ 

dpi 

' dz ' 


(4.3c) 


(4.3d) 


(4.4a) 

(4.46) 

(4.4c) 

(4.4d) 


Here the superscript * on a quantity denotes the complex conjugate. 

We note that the TS interaction appears at leading order in (4.3c, d). Bennett et al., 
(11), examined the problem with O(l) a- variations but we analyse the limit — > oo to 

obtain an asymptotic description of the long wave-short wave vortex-TS interaction. 

If /? «C 1, but of sufficient size that subsequent expansions in ascending powers of /? do 
not disrupt the leading order equations (4.3, 4.4) we find that on a z-lengthscale of 0(/?) 
the leading order vortex and TS quantities assume the forms 


«o — “o + • • • + cos ^ | + |/d 2 ip i + • • ■) cos > (4.5a) 

Vo = |/? _1 (itf + • • •) cos | + | O'? + • • •) cos (^jfj + • • • J , ( 4 - 56 ) 

two = | (Wq + . . .) sin J + {/? (w? + • . .) sin | j^J + . . .| , (4.5c) 

Po — P 4 Po + • • • + j/? 1 (Pq + • • •) cos | + j/? 2 {P\ + • • •) cos + * • • j , ( 4 -5d) 
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and 


ui = P~ 2 ut + . . . + j/T 1 ^ cos(i) +...}, (4.6a) 

«i = P~ 2 VJ + . . . + {/T 1 ^ cos ^ , (4.66) 

w i = pw} + . . . + sin + . . . j , (4.6c) 

Pi = P~ 2 p} + . . . + ^PPi cos + . . .| , (4.6d) 

where all the unknowns are function of y. 


We substitute these expansions into (4.3) & (4.4), together with the relevant expansion 
T = T 0 /? -4 + . . . for the Taylor number. On equating powers of /? associated with each 
Fourier coefficent we get the solution set 


rt 


, du o 
~dy' 


Uo = A^, V} = -iccAuo, Po = p- a 2 A f ( u 0 ) 2 dy , 

J o 


(4.7) 


where the displacement A and pressure term p are to be found. Further, we find that 

■M ?-¥$)■ (It-f )• <*■•> 


Additionally, the vortex terms Uq, Vq, Wq satisfy 


Vo' + ( 

and the mean flow term U 0 satisfies 


0 + 2 dy 0 ~ ’ 

2 ToSo- 16 a 2 |A| 2d “°^ 


dy 


Uq = 0 , 


d 2 u 0 


+ 


D — S5® W) *)- 


Now (4.96, c) are consistent only if 


(«•->£) 5 


= 1 , 


(4.9a) 

(4.96) 

(4.9c) 

(4.9d) 


(4.10) 


where we have put A = 8a 2 |A| 2 . 
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(4.11) 


The TS- vortex interaction requires that the mean flow satisfies (4.10) and then (4.9d) 
fixes the vortex. The TS wave is governed by the solutions (4.7, 4.8) and is fully deter- 
mined by considering the wall layers, see below. The consistency condition (4.10) is a 
generalisation of (3.6a), the additional term stemming from the presence of the TS wave. 
However, a major difference in the solutions of these two equations exists. Previously, for 
(3.6a), given u 0 at some point y* say, the solution is uniquely defined, whereas now we 
have 

dti o 7o§o ± \J (To«o) 2 — 4A 

dy 2A 

with a consequent ambiguity in the solution. This phenomenon will be found to play a 
crucial role; for the present we take the negative sign in (4.11) since when A — ► 0 this 
retrieves the analysis of §3. 

We may integrate (4.9d) once to obtain 

dti o _ 1 duo / « \ 2 

_ + = ___ (K 0 ‘) +*, 

where 6 is a constant and, as in the previous section, we expect the vortex activity to be 
concentrated in a region (7) lying between j/i and y 2 say. The condition that | V^ 1 1 
these positions together with (4.11) yields the relationship 


0 at 


To 
2A L 


u 0 - ^/(T 0 u 0 ) 2 - 4 A = D(y x - y 2 ). (4-12) 

L J v=vi 

We anticipate shear layers at t/x and y 2 as before, again of thickness 0(/?i) and which 


we examine now. 


§4.1. The shear layers. 

In the shear layer (I I a) centred on yi, we define the O(l) co-ordinate £ = /? _ ?(y — yi) 
and the expansions (4.5) & (4.6) are modified to give 


u 0 = uj + 0*u\ + /?uj + P*ul + . .. + j/3* (Uqi + P*Uu + . ..) cos 0}* 

= j/? - ® (koi + P J Vn + . ..) cos ^ | + . . . , 

= {/H(w 0 i + /?*w n + ...) Sin (j^)J + ..., 

= 0 4 Po + 0 3 Px + . . . + |/3 2 (Pol + 0* Pn + • • •) cos ^ | + 


Vo 


Wo 


(4.13a) 

(4.136) 

(4.13c) 

(4.13d) 
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and 


u, = r 2 ui +r*ul , + ...+ jr* (yl + fiul + . . .) «* (i) j + . . . , (4.140) 

vi = P-'VI + r *V„t + . .. + |r 1 (V,t + PWX + . ..) cos (i) } + .... (4.146) 

Wi = pwi + . . . + |r i (w^ + piw}, + . . .) sin (|) j + . . . , (4.14c) 

Pl = r !, pl+r‘P^ + --- + ip‘ (p 1 t „ + /j>p 1 t i + ...)cos^j| + ... . (4.i4rf) 

These forms follow from (4.5)-(4.11) and we substitute these expressions into (4.3) Sc 
(4.4). As in the previous section, it is found that the mean flow terms Uq,Ui satisfy 

d?ul _ d 2 u\ _ 

~d^~ = ~d? = ' 

A solvability condition for the fundamental vortex term V 0 i is derived by following the 
procedure given in §3 and determining circumstances under which the equations governing 
U\i and Vn admit a solution. The algebraic manipulation needed now is more complicated 
than before because of the presence of the TS wave. However, eventually we find that 

(3 - + [iw + - 1)] evil = ^(1 - AkM + fVm, (4.15) 


where Hi = ^ and / is some constant. Again, we have an unsealed form of the 

“V Vl 

second Painleve transcendent, and the solution has the desired property that Vox — ► 0 as 
£ — ► — oo and Voi = 0(£s) as £ — ► oo, as long as 

= *(1 - Arf) > °- (4 ' 16o) 

An entirely analogous analysis is appropriate for the shear layer ( lib ) positioned at y = 
y 2 and the resulting Painleve problem is precisely (4.15) but with replaced throughout 
by /x 2 = . Again, the fundamental vortex term V 0 i decays exponentially as we move 


away from region (I) if 


„ _ T„p.\ + DM - 1) „ „ 
o 2 — ; — ^ <- u. 


(4.166) 


/i 2 (l - Xnl) 

The conditions (4.16) may only be checked by numerical computations and so at this 

t t 

stage we assume that they hold. The fundamental terms U { 0 , . . . , 0 associated with the 

TS wave in expansions (4.14) are expressible in terms of the vortex coefficients Uqi, . . . , Pqi- 


14 



The exponential decay of these latter terms imply a similar behavior for the former and 
so outside regions (I) & (I la, 6), the implied velocity and pressure fields develop as 


A 

u 0 - u 0 + . . . , 

u 1 = p- 2 u% + ..., 

(4.17a, e) 

+ 

o 

II 

O 

+ 

to 

n 

1 

II 

H 

(4.176,/) 

too = o + . . . , 

Wi = PW$ + . . . , 

(4.17c, g) 

Po = P~% + • • • , 

Pi = P~ 2 Pq + . . . , 

(4.17 d,h) 


to leading order. We find that, as in regions ( IIIa,b ) examined in §3, the mean flow 

A 

quantity tto is simply given by 



(4.18a) 


and here the TS parts Uq,Vq 


satisfy 


Uq = A^p-, Vq = -iccAZ 0 . (4.186) 

ay 

Finally, to complete the flow structure which thus far has been identical to that given 
in §3, we add wall layers at y = 0, 1. In the layer at y = 0, we define the 0(1) co-ordinate 
by y = e 2 Y and the fluid velocity and pressure expand as 

u = e 2 U ol + . . . + ^e 6 /? 2 U i lE + . . . + c.c^j (l + 0(e 7 )) , 

v = (pP 2 V nE + . . . + c.c.j ^1 + 0(e 7 )J , 
w = (e 7 /3fr lL E + . . . + c.c.'j (l + 0(e 7 )) , 

p = —Dxe 7 + (fi~*PoL + ...) + (p/3~ 2 P 1L E + . . • + c.c.^ (l + 0(e 7 )) . 

Substituting these expressions into the governing equations (2.2) leads to the usual 
lower branch TS wave eigenproblem, see for example (8). We find that the TS pressure 
on the lower wall Y = 0 is given by 

p _ A\lAi{l L ) 


where Ax, is the shear of the mean flow u 0 at the wall, Ai is the Airy function (see (19)) 
and = — Mft/(aAi)?. Consideration of the upper wall leads to a similar TS pressure at 
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y = 1 and using equation (4.7c) for the leading order TS pressure across the whole of the 
channel, we finally obtain 


.'iaf ( /"(Bo) V) = 7 — 

Vo K 01 y ) ( fP° AilP\At 


^u-Ai (\u) 


<C A 'i 


(4.19) 


as the eigenrelation to determine the neutral wavenumber a and frequency fl. In this 
formula, \ v is the shear of the mean flow at the upper wall and \u = —i^Cl/(a.'\ u )l. 

We now have sufficient knowledge so as to completely specify the problem. In the core 
region the leading order mean flow satisfies (4.10) and (assuming that (4.16a, b) hold) has 
continuous first derivative at the shear layers y x and y 2 . These shear layer positions satisfy 
(4.12) and outside the core the mean flow is determined by (4.18a). The imposition of zero 
velocity conditions at the walls now (in principle) yields the complete mean flow across 
the channel to within an arbitrary constant. This degree of freedom is fixed by choosing 
the pressure gradient D so as to ensure constant non-dimensional mass flux and, finally, 
(4.19) yields the wavenumber a and frequency Cl of the TS wave. 

A A 

Analysis of the governing system reveals that it is convenient to define A by A = ATo. 
The problem is then easily formulated in terms of the parameter groupings TqD 2 and A. 
If we write the mean flow in region (Ilia), (where (4.18a) holds) as 

U 0 = D ($y - O , (4.20a) 


u 0 = D ^$y - yJ > (4.20a) 

where $ is an initial guess for the shear at the lower wall, then the position y x of the lower 
shear layer satisfies 

*V.-y-X(*-»i) = (4.206) 

It has already been noted that an ambiguity appears to exist as to which sign to take 
in the expression (4.11). However, it is easily shown that if the positive sign is taken, 
then it is not possible for (4.10) to hold for any y > y x ; consequently no upper shear layer 
could exist and the flow structure described would not be valid. We are thus led to the 
conclusion that the negative sign should be taken in (4.11), and if we write /2 0 = DU, the 
solution of (4.11) is of the form 


ToD’p TqD 2 
4 4 


U 2 - 


T 0 D 2 


ffoD 2 

' iA 

4A 


U + \ U 2 


4A 

TqD 2 


y + const., (4.20c) 


where the constant is chosen to ensure continuity in u 0 across y x . Equation (4.20c) together 
with (4.12) which now assumes the form 


1 jj _ fr 2 _ 4A 

2 A U \ U TqD 2 


= Vl- V2, 


(4.20c!) 
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yields y 2 and the mean flow and its derivative at y = in- 
equation 

d 2 u 0 


For y > y 2 , we revert to the 


dy 2 


= -A 


which is integrated to the upper wall, y = 1. A simple Newton iteration is applied to the 
guess for the lower wall shear value $ to ensure that §0 vanishes at y = 1. The constant 
D is then fixed so as to give constant mass flux and (4.19) gives a and Cl. 

The computations were initially performed by fixing T 0 D 2 and then repeatedly solving 

A A A 

the problem for increasing values of A from A = 0. We note that the limit A — » 0 
corresponds to the problem considered in §3 where no TS wave is present. As in the work 
presented there, no solutions were found for T 0 D 2 < 24a/ 3. For small values of T 0 D 2 up 

A. 

to approximately 80, as A was increased the shear layer positions yi and y 2 monotonically 

A 

increased and decreased respectively until a critical value of A was reached (dependent 
upon T 0 D 2 ) at which stage yi and y 2 merged. At this point, the vortex disappeared and 
the problem reduced to that of a channel flow with a small TS wave imposed with the 
resulting neutral wavenumber and frequency a = 5.7834, fl = 7.4016. This critical value 

A A 

of A, say A c is determined by demanding that equation (4.206) with $ = 1 has a double 

A A 

root at y = yi. For A > A c , the only solution which exists is plane Poiseuille flow with no 
vortex and a linear 2D-TS wave of arbitrary amplitude imposed. 

For values of T 0 D 2 >« 80, a different behavior was observed. As A was increased, 
but prior to yi and y 2 merging, the solution apparently vanished, although solutions could 
again be generated in a neighborhood of \ c (TqD 2 ) defined in the previous paragraph. It 
was noted that as this breakdown position was approached the value Si defined in (4.16a) 
approached zero and, in turn, this ensures that the two possible values of (^-) at y = y x , 
given in (4.11), merge. A detailed description of the results of these and the subsequent 
numerical solutions will be given in §5 where the significance of the results will also be 
discussed. 

If \ B i(T 0 D 2 ) and X B2 (TqD 2 ) are respectively the values of A at which breakdown occurs 

A A 

as A is increased from zero and decreased from A c , at these values the Painleve problem 
(4.15) at y = yi clearly becomes inapplicable. The shear layer is again of thickness 0(/?a) 
but the fluid quantities now behave as 


f it t 
Uq — Uq + P Zu l + P u 2 + • • • + 


Vo 


Wo 


{/?* {Uoi + (3>Un + ...) cos 
= |r l (Voi + ^Vn + ...) cos + 

= |r ! * (Woi + 0*W U + . . •) sin ) + . . . , 


+ 


(4.21a) 

(4.216) 

(4.21c) 
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Po — P 4 Po + P 3 Pi+-**+|^ 6 (P 01 + P>Pu + ... 

) cos (|)} + --' 

(4.21d) 

«1 = p- 2 ul + p-iui + . . . + |r 1 (ut + piul + ... 

) cos (i)} + -’ 

(4.22a) 

v\ = P 2 V 0 o + P 3 Vqi + . . . + |/? * (Vjo + P 3 V n + . . . 

) c °s(i) } + ..., 

(4.226) 

ttfi = pwl + . . . + j/H (wtfo + P*w} x + . . .) sin 

6)} + -- 

(4.22c) 

Pi = P~ 2 pI + P^Pox + • . ■ + {/?* (Pio + P*Pn + ■■■] 

) cos (^)| + > "’ 

(4.22 d) 


where these expansions follow from the solutions in region (/) and the breakdown assump- 
tion 

Ton* = -^(1 — ^Pi)> 


from (4.16a). We may retrieve a solvability condition for the fundamental vortex term Voi 
which takes the form 

Vn 


*Vo. To 

~W + 2A fVo1 


'01 

128 


(ySt-f.) , 


(4.23) 


where /* is a constant, c.f. (4.15). This equation has a solution with Voi =0(^1) as 
£ — > oo (which matches with the behavior of the quantity Vq 1 in the core as y — *• t/i) 
and exponential decay as £ — ► — oo, and so the overall flow structure elucidated already 

AAA 

remains valid. We note that at A = Xbu^bi the structure of the upper shear layer at y 2 
is not affected (i.e. < 0 in (4.166)) and so the modified version of the second Painleve 

transcendent is still the governing equation there. 


A A A 

§4.2. A model problem for the case A B1 < A < A B2 . 


We conjecture that a new solution structure must govern the flow for the range of A 
between the breakdown points. To examine this case, we considered the simpler model 
problem given by 



(4.24a) 



y- 



(4.246) 
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with 0 < t <1, u(0) = u(l) = U(0) = 27 (l) = 0 and where A is some constant to be 
determined. The relevance of this system to the problem at hand is thus; here we expect 
to have a central region where U = 0(1) and 


1 * ( du\‘ 

T~ t + * [d^J 


du 

~ U T = °» 

dy 


(4.25) 


(c.f. (4.10)), bounded by thin regions in which U changes rapidly to zero (c.f. shear regions 
{II)) and in the remainder of the channel ^ » A — y, (c.f. regions ( IIIa,b )). Numerical 
solution of (4.24) for e small but non-zero revealed the existence of a non-trivial solution 
for all 0 < A* until the critical value of this parameter was reached at which the ‘shear 
layer’ positions merged. This further suggests that a new structure for the asymptotic 
problem comes into play at some stage. 

For asymptotically small e, the structure proposed in §§4, 4.1 is valid for this model 
problem for 0 < A* < \* B i{T*) and ^ sa y where \ Bl [T*) and ^*B 2 i T *) are the 
breakdown points described before. For < A, < A^ 2 (T*), we suppose that a shear 

layer exists at y = yj, and that it has thickness 0(c), i.e., much thinner than the previous 
0(c?) size. If in the vicinity of y{ we define the 0(1) co-ordinate ip by 

y -y*\ = (4.26a) 

then the functions U, u behave as 

U = U 0 + . . . , u = <p 0 + eu 0 + , (4.26 b) 

locally. Here <p 0 = Ay{ — (yj 2 / 2) is the value of « at y = yj, and we further suppose that 
just below yj we have = <pi = A — yj and U — * 0. Putting (4.26) into (4.24) yields the 
system 

(4.27a) 


d?U o 

dip 2 


= U 0 



(4.276) 


Equation (4.276) may be integrated once to obtain the solution 

O’ ■ [I (rife) -«"■('*«)] 


which satisfies — ► 0 as U 0 — ► 0. 


These equations admit a solution in which ^ — *• <p\ (< <pi) as ip — ► oo (and so 
Uq — > const, from (4.27a)) where 

‘"(It) 


— + \*<pi<p{ — <po<Pi<Pi 


{<f>i ~ <f > 0 


= 0 . 


(4.28) 
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(4.29) 


Further, we fix <f>{ to satisfy (4.25) so that 


- + a. (ft) 1 - ow; = o, 

1 * 

and across the central region u is governed by the solution of 


du _ 1 

dy 2A* 





2 

U 2 - — - 

T* 


(4.30) 


subject to u{y{) = 4> 0 . 
(4.24 b)) by 


Finally, the upper shear layer position y% is determined (using 


y{ - y\ = 


2A* 


A 



(4.31) 


and for all the problems encountered in this paper, the shear layer at y 2 was found to be 
of the thicker 0(c?) variety with the usual Painleve-type governing equation. 

To summarize this new structure, we have shown that now a thinner shear layer may 
exist across which u is continuous and ^ suffers a discontinuity. In this layer, the ‘vortex’ 
term U is of size 0(1) in contrast to its much smaller magnitude in the case of the O(es) 
layer structure of §4.1. This term tends to zero as we move away from the core and to an 
0(1) constant as we approach the core. The calculation for A^ 1 (T*) < A, < A g 2 (T*) was 
tackled in the following manner. 

As before, the shear for the ‘mean flow’ quantity u at the lower wall ( A in (4.246)) 
was guessed. Previously, for the thicker shear layer, the position of the lower layer was 
determined by demanding that (4.25) be satisfied, but now this is no longer the required 
condition. Instead, we now solve (4.28) & (4.29) for <f>\ and <p J and hence determine y{. 
In turn, (4.30) & (4.31) are used to find yj, above which U — ► 0 exponentially and u(l) 
found using (4.24). A Newton iteration on A is then used to ensure that u(l) = 0. This 
calculation demonstrated that the solution structure proposed was the required one for 
*ii(r.) < A* < A^ 2 (T*), so that depending on the values of A* and T* one of the two 
distinct structures provided the solution to our problem. 

The result of this investigation of the model problem suggested that a similar phe- 
nomenon of a change in the asymptotic solution structure may occur in the vortex-TS 

AAA 

interaction problem under consideration. If for A^i < A < Xb 2 we seek a structure with 
the lower shear layer at yi of size 0(/?), then the solutions for the unknowns of (4.5) & (4.6) 
in the region of vortex activity suggest that in the shear region where now f = (y - 
we have , / \\ 

u 0 = ul + + . . . + | P ( Uqi + (dU\\ + . . .) cos ( — J | + . • • , (4.32a) 


20 



Vo 


Wo 


— j/3 1 (Voi + PV\\ + . . .) cos ^ | + . . . , 
= { rl 


{W 01 + PW 11 + —) sin (- 

P / 


+ . . . , 


Po = /? 4 Po + 


(' 


+ \ fi 3 (Poi + &P\l + • • •) COS 


and 


(f)} + - 

«i = r' 'tr & + . . . + |r 1 (trio + pat + . . .) cos f | + . . . 

«i = /r’Vot + . . . + jr* (v,t + «! + . . .) cos (|) \ + 

®1 = + • • • + {r 1 (^4, + /Wi + . . .) sin f i j 1 + 

pi = r ! 4,+ ••• + {(4,+ ppI , + ...) cos j + . . 


(4.326) 

(4.32c) 

(4.32d) 

(4.33a) 

(4.336) 

(4.33c) 

(4.33d) 


instead of (4.13) Sc (4.14). After substituting these expansions into equations (4.3) Sc (4.4) 
and some manipulation, we find that uj = const. (= (f>o say, in the notation of the model 
problem discussion) and that we have the solvability problem 


t 


d 4 V 0 i d 2 V 0i , s du\ 

-jgr ~ 2 ~d£T + Voi + TotoUoi = A(J7oi + V 0l -j± 


Voi~ 


d 2 V 


01 


d? 


— Uqi + 


d 2 U, 


01 


d? 


). 


duj A. i_ ^Voi 
Of = * ,+ 4 ’ 

d’Vo, l v (duf 

-5j5- - = -V 01 — 


(4.34a) 

(4.346) 

(4.34c) 

(4.34d) 


rt 


where = iaAV 0 i and <f>i is a constant. To meet the boundary conditions required 
we need to specify that the vortex decays exponentially below t/i so that Uqi, Vo„V? 0 — » 

0, — ♦ d>! as £ — > — oo. We additionally need and the vortex terms to approach 

non-zero constants as £ — ► oo. That (4.34) admits such a solution is not immediately 
obvious, but this system bears close analogy with the model problem already discussed. 
If the full problem is derived with a small parameter multiplying the pressure in the axial 
momentum equation then, if this parameter is allowed to tend to zero, we recover the model 
equation. The structure of (4.34), together with the similarity with the model (4.24) and 
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the subsequent solution of that system suggest strongly that a solution of (4.34) with the 
required behavior exists. 

If the mean flow quantity uo inside the activity region is defined by &o = DU as before, 
then here 

dU_ _ J_ 
dy 2A 

and y 2 is given by 


A 

U -\ U '~ 3 ^ 5 ? - *>• < 4 - 356 ) 

J y=y 2 

where (4.28) &: (4.29) are to be satisfied with T* replaced by T 0 D 2 in these equations. 

This modified structure was then used to numerically solve the interaction problem for 
Abi < A < Ab 2 - The results and some discussion are now given in §5. 

§5. Results and discussion. 

We now discuss the results for the strongly nonlinear TS-vortex interaction which can 
be found using the scheme outlined in the previous section. We note that the numerical 
scheme used allows for the possibility that the lower shear layer bounding the region of 
vortex activity to be either of thickness 0(/?s) or 0(0) depending on the properties of the 
core equation (4.10) at yi. 

In figure (4) , we have shown the dependence of D on the Taylor number To for several 
values of A = A/To which characterizes the size of the TS wave. We note that all the 
curves terminate on the line D = 2 at which point there is no vortex activity and the flow 
reduces to plane Poiseuille flow with a superimposed two-dimensional TS wave of arbitrary 
amplitude. The wavenumber a and the frequency fl of course vary on each of the curves 
and on D = 2 a and fl have the linear neutral values for plane Poiseuille flow, namely 
a = 5.7834, fl = 7.4016. In figure (5), we have shown the development of two of these 
curves at higher values of To- The shapes of these curves are typical of the results we found 
for all A, thus when T 0 — * o o, we have D — ► oo but D/T 0 — * 0. 

In figure (6), we show the variations of the shear layer positions yi and y 2 with To for 
the cases of A = 0.01 and A = 0.2. The values of yx which correspond to the Painleve 
transition layer as opposed to the structure discussed in §4.2 are indicated by the part of 
the curve which lies between the crosses marked on figure (6). The values of yi and y 2 
coalesce when D — » 2 in which case the basic state in the absence of the vortices is simply 
plane Poiseuille flow. The mean velocity profiles corresponding to figures (5u), (6tt) are 
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shown in figure (7) for a few representative values of To. For A = 0.2 vortices may only exist 
for T 0 > 19.06 and the Painleve shear problem at y x switches to the other structure with 
a much thinner shear layer when To « 21. In figure (7*), for which T 0 = 19.98 the mean 
flow has a continuous derivative at yi and the overall profile is not much changed from the 
plane Poiseuille case T 0 < 19.06. In figure (7tt), To = 23.52 and the discontinuity in the 
derivative of the mean flow at yi is clear now that the very thin shear layer structure has 
taken over. For even higher Taylor numbers (To = 26.61 in fig. (7 tit)) this discontinuity 
becomes increasingly more marked. 

It should be stressed that figures (5) & (6) do not describe directly the bifurcation 
structure of a finite amplitude TS wave where T 0 or D varies. In an experimental setting, 
the TS wave is often generated by a vibrating ribbon so that as it moves downstream its 
frequency is constant. The bifurcation picture for constant frequency TS waves can be 
deduced from figure (4) if the lines of constant fl are plotted. However, it is easier to 
calculate the bifurcation curves by modifying the algorithm described in §4 to allow for an 
outer loop in which we can iterate on T 0 D 2 to make fl constant. 

In figure (8), we show the bifurcation pictures for five fixed frequency TS waves. We 
take a|A| as a representative TS amplitude and show how this quantity varies with T 0 for 
the frequencies chosen. We see that in each case there is a supercritical bifurcation to a 
finite amplitude mixed vortex-TS state. At the bifurcation point on the horizontal axis the 
TS wave has zero amplitude and the basic state is the fully nonlinear vortex flow described 
in section 3. 

Solutions for the nonlinear vortex state have been shown to exist only for Taylor 
numbers To < 32. Calculations performed for increasingly large values of fl (curves for 
fl = 1200 & 1600 are shown in figure (8)) suggest that for Taylor numbers To > 32 TS 
waves of arbitrarily small amplitude are possible. As fl — > 7.4016 there is a limiting 
bifurcation curve (dotted in figure (8) and labelled C{) corresponding to the linear neutral 
frequency for TS waves in Poiseuille flow. We note that it is straightforward to show that 
for large To the amplitude of the TS wave on C\ is proportional to T 0 . All the bifurcation 
curves lie to the right of C\ and joins the horizontal axis at a point 6\/3 < T 0 < 32. 

Now let us turn to the physical implications of the above results. Firstly suppose that 
the Taylor number is fixed. If To is less than the linear critical value 6\/3, then the only 
possible state is plane Poiseuille flow. If To lies in the range 6\/3 < To < 32, then there 
exists either a mixed TS-vortex state, a finite amplitude vortex state or plane Poiseuille 
flow. Since the mixed state bifurcates supercritically from the vortex state, it is presumably 
the only available stable solution of the vortex-TS interaction equations. In this Taylor 
number range, there is an upper limit to the size of the TS wave which the vortex flow can 
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tolerate. Thus, when C\ is approached the vortices decay and the azimuthal velocity field 
reduces to plane Poiseuille flow. Where C\ is crossed the vortex-TS interaction equations 
no longer determine the size of the TS wave which will then be fixed by a weakly nonlinear 
theory based on the Stuart-Watson method. Thus, in the Taylor number region where a 
strongly nonlinear vortex flow exists, the imposition of a TS wave of gradually increasing 
size ultimately destroys the vortex state and sets up a flow comprising plane Poiseuille 
flow and a two-dimensional TS wave. At Taylor numbers beyond 32, the range of values of 
a|A| where a mixed vortex-TS flow can exist lies below the curve C\ of figure (8). Outside 
this region only plane Poiseuille flow is possible. Hence at Taylor numbers greater than 32, 
finite amplitude vortices can exist only in the presence of TS waves. Thus, in this regime 
the TS waves are responsible for the existence of any vortices in the flow. 

Our results therefore suggest that experimentally there can be a significant difference in 
the response of the flow to a TS wave dependent on the wall curvature. At relatively small 
curvatures, small amplitude TS waves will equilibrate to a size fixed by its frequency. If that 
frequency is decreased, then the TS wave size increases until the vortex flow is destroyed. 
At higher curvatures no vortex flow (with the given prescribed wavelength presumably 
introduced by wall roughness or turbulence in an experiment) can exist unless a TS wave 
is also present. If the TS amplitude is decreased, the vortices are again destroyed. Hence, 
at any Taylor number a sufficiently large TS wave introduced into the flow will destroy 
any vortex flow. 

It is of interest to note that the calculations suggest that no mixed vortex-TS state 
can persist at Taylor numbers below the linear critical value of 6\/3. The work of Hall & 
Smith (13,14) suggests that the vortex-TS interactions at larger vortex wavelengths can, 
even without curvature, lead to singularities developing simultaneously in both time and 
the spanwise co-ordinate. Our results suggest that the transfer of energy is possible on 
smaller scales and so the results of Hall & Smith do not continue indefinitely. 
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FIGURE CAPTIONS 


Fig. (1). Positions of the lower and upper shear layers t/i, yi as a function of the Taylor 
number To- 

Fig. (2). Dependence of the Taylor number T 0 on the pressure gradient D when the 
constraint of constant non-dimensional mass flux is imposed. 

Fig. (3). Graph to show the dimensional mass flow rate as a function of the pressure 
drop in a realistic experiment. In the absence of vortices the relationship is linear (solid 
line). When vortices are present, one of the broken lines is followed (which of the curves is 
determined by the physical parameters of the apparatus) which show that as the mass flux 
slowly rises an increasingly large change in pressure gradient is required. Indeed, for any 
particular case the possible values for mass flow rate are bounded above and these bounds 
are indicated for the lower two curves in the figure. 

Fig. (4). Graph of the pressure gradient D against Taylor number T 0 for a variety of 

^ A 

values of A = 8a 2 |A| 2 /T 0 . The lines (l)-(10) each follow a path of constant A up to the 
point where T 0 D 2 = 1000, and are given by (1) A=10~ 4 , (2) A= 10 -2 , (3) 5 x 10 -2 , (4) 0.1, 
(5) 0.2, (6) 0.3, (7) 0.4, (8) 0.5, (9) 0.8 & (10) 1.0. The dotted line D = 2 corresponds to 
the problem where no vortices are present and the flow is plane Poiseuille with a linear TS 
wave of arbitrary amplitude imposed. 

Fig. (5). The development of two of the curves of fig. (4) for larger values of T 0 . (i) 
A = 0.01 (curve (2) in fig. (4)), (ii) A = 0.2 (curve (5) in fig. (4)). 

Fig. (6). Position of the shear layers yi and t /2 as functions of the Taylor number To for 
the cases (i) A = 0.01, (ii) A = 0.2. The positions of t/j which correspond to the Painleve 
shear layer problem lie between the crosses marked on the lower curves. 

Fig. (7). Mean flow velocity profiles for A = 0.2 when (i) T 0 = 19.98, (ii) T 0 = 23.52 & 
To = 26.61. The dotted lines represent the positions of the shear layers. 

Fig. (8). Constant frequency curves in Taylor number/amplitude space. The values on 
each line are (i) Cl = 10, (ii) Cl = 25, (iii) 0 = 36.5, (iv) Cl = 1200 k (v) fi = 1600. The 
dotted line (C7 X ) corresponds to the limiting case Cl — > 7.4016. 
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